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Abstract 

Oscillons are extremely long-lived, oscillatory, spatially localized field configurations 
that arise from generic initial conditions in a large number of nonlinear field theories. 
With an eye towards their cosmological implications, we investigate their properties 
in an expanding universe. We (1) provide an analytic solution for one-dimensional 
oscillons (for the models under consideration) and discuss their generalization to three 
dimensions, (2) discuss their stability against long wavelength perturbations, and (3) 
estimate the effects of expansion on their shapes and lifetimes. In particular, we discuss 
a new, extended class of oscillons with surprisingly flat tops. We show that these flat- 
topped oscillons are more robust against collapse instabilities in (3+1) dimensions than 
their usual counterparts. Unlike the solutions found in the small amplitude analysis, 
the width of these configurations is a nonmonotonic function of their amplitudes. 
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1 Introduction 



A number of physical phenomenon from water waves travehng in narrow canals [T] , to phase 
transitions in the early Universe [2] exhibit the formation of localized, energy density config- 
urations, even without gravitational interactions. The reason for their longevity are varied. 
Some configurations are stable due to conservation of topological or nontopological charges, 
while some are due to a dynamical balance between the nonlinearities and dissipative forces. 

Relativistic, scalar field theories (with nonlinear potentials) form simple yet interesting 
candidates for studying such phenomenon. Some well-studied examples include topological 
solitons in the 1 + 1-dimensional Sine-Gordon model and nontopological solitons such as Q- 
balls [3j . The Sine-Gordon soliton is stationary in time whereas the Q-balls are oscillatory in 
nature. Both have conserved charges which make them stable (at least without coupling to 
gravity). This paper deals with another interesting example of such localized configurations 
called oscillons (also called breathers). Like the Sine-Gordon soliton, they can exist in real 
scalar fields, and like the Q-balls they are oscillatory in nature. Unlike both of the above 
examples they do not have any known conserved charges (however, see 0] for an adiabatic 
invariant). In general they decay, however their lifetimes are significantly longer than any 
natural time scales present in the lagrangian. Along with their longevity, another fascinating 
aspect of oscillons is that they emerge naturally from relatively arbitrary initial conditions. 

Not all scalar field theories support oscillons. In the next section we discuss the require- 
ments for the potential. Here, we note that the requirement is satisfied by a large number of 
physically well-motivated examples. For example, the potential for the axion, as well as al- 
most any potential near a vacuum expectation value related to symmetry breaking, support 
oscillons. Oscillons have also been found in the restricted standard model SU{2) x U{1), 

Oscillons first made their appearance in the literature in the 1970s |H]. They were sub- 
sequently rediscovered in the 1990s [9j. Oscillons are not exact solutions and (very slowly) 
radiate their energy away. The amplitude of the outgoing radiation (in the small amplitude 
expansion) has been calculated by a number of authors, see for example [IHl [TTl |12] . Charac- 
terization of their lifetimes and related properties using the "Gaussian" ansatz for the spatial 
profile was done in |I3J (also see references therein). The importance of the dimensionality 
of space for these objects has been discussed in [HI [T^ . 

Their possible applications in early Universe physics has not gone unnoticed. For ex- 
ample, they could be relevant for axion dynamics near the QCD phase transition [16]. The 
properties of oscillons in a 1 + 1-dimensional expanding universe (in the small amplitude 
limit) have been discussed in [T71 [7]. Their importance during bubble collisions and phase 
transitions have been discussed in p!B]. In [TH], interactions of oscillons with each other and 
with domain walls were studied in 2-1-1 dimensions. 

In this paper, we point out what is required of scalar field potentials to support oscillons. 
We then derive the frequency as well as the spatial profile of the oscillons for a class of 
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models under consideration. We show that the spatial profile can be very different from a 
Gaussian, an ansatz often made in the literature. In particular, we derive the nonmonotonic 
relationship between the height and the width of the oscillons, and discuss the importance 
of this feature for the stability of oscillons (see [20] for a somewhat related analysis for Q- 
balls). To the best of our knowledge, this has not been done previously in the literature in 
the context of oscillons. We consider the stability of oscillons against small perturbations, 
mainly with spatial variations comparable to the width of the oscillons. We also comment 
on a possible, narrow band instability at higher wave numbers. 

Oscillons could have important applications in cosmology, especially in the early Universe. 
With this in mind, we discuss the changes in the profile and the loss of energy from these 
oscillons due to expansion. 

The properties of oscillons can depend significantly on the number of spatial dimensions. 
In this paper, for simplicity we always start with 1 + 1-dimensional scenarios where analytic 
treatment is often possible. We then extend our results to the physically more interesting 
case of 3 + 1 dimensions, analytically where possible and numerically otherwise. We extend 
previous analysis to the interesting "flat-top" oscillons because of our new systematic method 
for capturing the entire range of possible amplitudes, while still using the methods from the 
small amplitude expansion. Our expansion, which can also be thought of as a single frequency 
approximation, allows us to use results existing in the literature for time periodic, localized 
solutions. 

Although interesting as classical solutions, a quantum treatment can lead to changes 
in oscillon lifetimes [21J. Another question worth investigating is the stability of oscillons 
coupled to other fields. These two questions are beyond the scope of this paper. 

The rest of the paper is organized as follows: In Sec. |2]we give a brief overview of oscillons 
in scalar field theories. In Sec. [3] we introduce a simple model that is used throughout the 
paper. Section |4] deals with the derivation of the shape and frequency of the oscillons in the 
absence of expansion. Section [5] focuses on linear stability of oscillons. Section [6] discusses 
the effects of expansion. Our conclusions and future directions are presented in Sec. [7j 

2 A gentle introduction to oscillons in scalar fields 

Oscillons are extremely long-lived, oscillatory, spatially localized field configurations that 
exist in a large number of nonlinear scalar field theories. We find it convenient to visualize 
an oscillon as a spatially localized, smooth envelope of the field value oscillating with a 
constant frequency. To get a heuristic understanding of what kind of the potentials support 
oscillons, let us consider the equation of motion for a 1 + 1-dimensional scalar field: 



d 



Dcp - V'icf) = 0, 
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where V'{(f) — )■ m'^Lp as 0. Let us approximate the oscillon field configuration as 

(p{t,x) ~ $ (a;) cos [cut]. For a locahzed configuration, as we move far enough away from the 
center (whereby the nonhnearity in the potential is irrelevant), we get 

_ + _ 0. (2) 

Again, because we are looking for a smooth, localized configuration, we must have u"^ < m^. 
For regions near the center of the configuration, we expect 9^$ < for the lowest energy 
solutions. This means that for the equation to be satisfied, 

V'{^)~m^^ < 0. 

Thus for oscillons to exist in potentials with a quadratic minimum, we require V'{ip) < m?ip 
for some range of the field value. 

It is not too difficult to think of physically motivated potentials satisfying this require- 
ment. For example, the potential for the QCD axions V{ip) = m?p [1 — cos(<y9//)] where / 
is the Peccei-Quinn scale and m is the mass, or any symmetry breaking potential expanded 
about it's vacuum expectation value. Both potentials "open up" a little when we move away 
from the minimum. 

The above (heuristic) argument does not provide a reason for the longevity of oscillons. 
For oscillons, their shape, which determines their Fourier content, guarantees that the am- 
phtude at the wave number of the outgoing radiation is exponentially suppressed (at least 
for the small amplitude oscillons). For details see [21] and the subsection on radiation in 
this paper. 



3 The model 

We begin with the action for a real scalar field in a ci + 1-dimensional, spatially fiat, homo- 
geneous, expanding universe {h = c = 1): 



Sd+i = J {adxYdt 



where 



(3) 



V{^) = -mV--^v'+l^\ (4) 

a{t) is the dimensionless scale factor and A, (? > 0. As we have discussed, what is crucial for 
the existence of oscillons is 

\/'((^) -mV < 0, 

for some range of the field. The potential [Eq. (|4])] is the simplest model which captures 
the effect we wish to explore, namely, a nonmonotonic relationship between the height and 
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width and its implication for stability. However, apart from detailed expressions, our results 
are general and not restricted to the particular shape of the potential. 

We find it convenient to work with dimensionless space-time variables as well as fields. 
Using {t,x) — 7- m'^{t,x) and (p — t- m\~^^'^Lp as well as (7 — ?■ {\/m)'^g the action becomes 



Sd-^ 



m 



3-d 1 



/ {adxYdt 



2a? 



(5) 



with 



J- 2 J- 4 , 5- 6 
2^ -4^ +6^- 



The classical equations of motion are given by 



-if + Hdt^p + 



0, 



(6) 



where g is the only free parameter in the potential and H = a/ a. We will be concentrate 
on the case where (7^1. This gives a controlled expansion in powers of g'^^"^^ which allows 
us to derive an analytic form for the profile. Our approach is similar to the small amplitude 
expansion, with the important difference that it captures the entire range of amplitudes for 
which oscillons exist. Moreover, in our analysis we show that the flat-top oscillons are stable 
against small amplitude, long wavelength perturbations on time scales of order g. 



4 Oscillon profile and frequency 

In this section we derive the spatial profile of the oscillons in our model in a 1 + 1- and 
3 + 1-dimensional Minkowski universe. We include the effects of expansion in Sec. [6j For 
simplicity, we begin with the 1 -f- 1-dimensional case. 

4.1 Profile and frequency in 1 + 1 dimensions 

The equation of motion is 

a,V-5.V + ^-^' + ^?^' = 0. (7) 
To extract the oscillon profile, we introduce the following change of variables: 

^{t,x) = -^(f){T,y) , 

(8) 

y = x/y/g, 



4 



where 



Here, characterizes the change in frequency due to the nonhear potential. We define 
$0 = 0(0,0) = y/g^o and choose dt^{^,x) = dr4>{0,x) = 0. Note that a and $o are not 
independent of each other. Their relationship will be determined from the requirement that 
the solution is periodic in time, smooth at the origin and vanishing at spatial infinity. With 
the change of variables (Isj), and collecting powers of g, the equations become 



+ <P^] = 0[g-'/% 



9^0 + + g-' {-a'dlc^ - - 
Let us consider solutions of the form 

v) = y) + y) + --- 

Again collecting powers of g, we get 

0, 



d% + 01 

a^03 + 03 



a2920i + 9j0i + 0? - 0^ 



The first equation in (11) has a solution of the form 

0i(r,|/) = ^{y) COST. 



(9) 



(10) 



(11) 



(12) 



To determine the profile $(?/), we look at the second equation in (11). Substituting 0i(r, y) 
into this equation, we get 



a^03 + 03 



^48 



COST, 



+ 



1$3 _ A$5 

4 16 



COS 3r 



1 

16 



(13) 



$5 cos5r. 



We are looking for solutions that are periodic in r. The term [. . .] cost will lead to a term 
linearly growing with r. Hence, we must have 



y 4 8 

This equation has a first integral, the "conserved energy" 

1 



0. 



(14) 



(15) 



where U ($) = — la;^^^ + — If we demand spatially localized solutions, we require 

Ey = 0. Furthermore, requiring that the profile be smooth at the origin, we must have 
(?y$(0) = 0. This immediately yields (also see Fig. [T]) 
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Figure 1: The above figure shows «^ which characterizes the change in frequency of oscillation 
due to the nonlinearities in the potential. The critical = a/27/160 can be obtained from 
the requirement that the nodeless solution is smooth and localized in space. Note that in 
1 + 1 dimensions, a is monotonic in $o- This is not the case in 3 + 1 dimensions. Frequency 
is measured in units of the m. 



= - (16) 

Note that there is a critical value etc = ^^27/160 (or $c = \/9/10) beyond which localized 
solutions do not exist. The profile equation becomes 

Integrating the above equation yields 



$(y) = $o./-^ ^-^^f;^, (18) 

^' 1 + 1* cosh [2Q;y]' ^ ^ 

where 

u = y/l - {a/acf, ^^^^ 
$0 = $cVl - u. 

We have introduced the variable < m < 1 which simplifies the appearance of the equations 
and controls the shape of the oscillons through a. We will come back to a more detailed 
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analysis of this solution, but first we solve for the second order correction to this solution, 



1$3 _ A$5 

4 16 



The solution with 03(O,y) = dr4>3{0,y) = is 



COS 3r 



1 

16 



-$ cos5r. 



(20) 



— (3$^ -4$^) COST + — (5$^ -4$3) cosSr + — $^cos5r. (21) 
96 ^ ^ 128 ^ ^ 384 ^ ^ 



The full solution becomes 

$3 

24g 



(r, y) = ^ COST + 



1 3 1 ' 
- (3 - 4$^) cos r (4 - 5$^) cos 3r H $^ cos 5r 



(22) 



Note that the corrections to 0i are strongly supressed for g ^ 1. Even for moderately large 
(7 ~ 5, the factor in the denominator is ~ 100, making (pi a rather good approximation. 
From now on we will mainly concern ourselves with 0i. 



Reverting back to the original variables (|8]), the solution for a < ac (equivalently, $0 < 
is 



l+U 



u cosh[2ax/y^] 



cos(a;t) + 



(23) 



where 



u 



$0 $ 



y/9 



(24) 



Here, ipo is the amplitude of the profile at the origin and scales as 1/y^. 



Let us now investigate the solution for the profile. Figure [2] (top left) shows this solution 
for different valued of a. Notice that as a approaches ac (equivalently, $0 u ^ Q), 

the oscillon profile begin to deviate from the "sech" profile and has a flat top. Given this 
solution, one can derive the width of the oscillon as a function of its height. Defining the 
width to be the x value where the profile falls by 1/e of its maximum 



1 2 
V^o a/3 



(1 + m)-i/2 pQgi^-i 



6^1 + u) 



u 



(25) 



As u — 7- we simply have Xe ~ l/v^o? which is consistent with the small amplitude analysis 
(see Fig. |3|. Meanwhile, u — )■ 1 yields a spatially uniform solution. 



We end this subsection by writing down an expression for the energy of these oscillons 

4 



En 



tanh 



3 1 



u 



1 - u 



1 



u 



+ 0[g 



-3/21 



(26) 
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Figure 2: The above figure shows the spatial profiles of oscillons for different values of the 
amplitude at the center. For $0 <^ $c = a/9/10 we get the usual sech-like profile, which 
is consistent with the small amplitude analysis. As $0 approaches $c, the oscillons become 
wider with surprisingly fiat tops. Unlike the 1 + 1-dimensional case, in 3 + 1 dimensions, we 
approach the fiat-top profiles from above. Distances are measured in units of the m~^. 



Note that as a — )■ (m — )• 1) we have E 

osc ~ 2^^2/3^90 whereas for — y ct^ — ^ 0), 
Eosc ^ 00. In the next subsection we extend our results to 3 + 1 dimensions. 



4.2 Profile and frequency in 3 + 1 dimensions 

In this section we extend the results of the previous section to a 3 + 1-dimensional Minkowski 
space-time. Although we are unable to obtain an analytic form for the profile, the important 
quahtative (and some quantitative) aspects of the solutions can still be understood. In 
particular, we derive a critical amplitude and frequency for which the solution becomes 
spatially homogeneous and argue that the relationship between the height and width is 
nonmonotonic. 

The equation of motion (assuming spherical symmetry) is given by 

- d^if - \(f + (f - (f^ + g^^ = 0. (27) 

We can follow the same procedure used in the previous subsection to arrive at the equation 
for the profile 

+ ^ap$ - + ^<l>^ - ^<l>^ = 0, (28) 

where p = r/^. This is where we first encounter the difficulty associated with three 
dimensions. We can no longer obtain a first integral due to the 2/p(9p$) term. However, we 
can still get a bound on a by requiring that the solutions are spatially localized (see [22] for 



8 




Figure 3: The above figure shows the nonmonotonic relationship between the width and 
height of oscillons in 1 + 1 and 3 + 1 dimensions. Note that as $o approaches $c = \/9/10, 
the oscillons become wider with flat tops. Unlike the 1 + 1-dimensional case, in 3 + 1 
dimensions, we obtain flat-top proflles when $o approaches $c from above. Distances are 
measured in units of the m,~^. 



an analysis of a similar proflle equation in the context of Q-balls). It is convenient to deflne 
an energy Ep, which in the absence of the {2/ p)dp<^ term, is a constant of motion: 

Ep=\{dp^f+Um, (29) 

where = — |q;^$^ + — With this deflnition the equation of motion takes on 

an intuitive form 

f^ = --(W- (30) 

dp p 

This means that as we move away from p = 0, we move from a higher Ep trajectory to a lower 
one. With the requirement that the solution is "localized" (more speciflcally, $ oc p~^e~°'P 
as p — > oo), we need £'p — )■ as p — )■ oo. Requiring that the solution is smooth at p = 
requires (9p$ = at p = 0. This implies that for a localized solution we must have Ep > 0. 
Equivalently, ?7($o) ^ 0, which in turn implies that a < = a/27/160. For this critical 
value we get a special solution [with $c(p) = v^9/10], which is homogeneous in space. 
For < Q < a^. wc get nonzero spatial derivatives. 

For each a in the range < a < etc, only special, discrete values of $o("') corresponding 
to the number of nodesjwill yield solutions that satisfy our requirement $ oc p^^e^"^ as 
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p — )■ oo. From these, the n = ones are the oscillon profiles we are looking for. The 
numerically obtained profiles are shown on the right in Fig. [2j 

From figure [3} it is easy to see that the relationship between the heights and widths of 
the oscillons is nonmonotonic. We know that for a <^ Uc (ie. $o ^ ^c), the usual small 
amplitude expansion yields solutions that have the property that their widths decrease with 
increasing amplitude. We also know that for a = ($o = ^c) the width will be infinite. 
Thus, as in the 1 + 1-dimensional scenario, we expect the width to be a nonmonotonic 
function of the central amplitude. This is indeed what is seen from the numerical solutions 
of the profile equation as shown in Fig. [3j Note that the width is a multivalued function of 
the amplitude beyond $o = ^c- Nevertheless, it still approaches the homogeneous solution 
via the fiat-top profiles. The multivalued relationship between a and $0 is shown in Fig. [T| 

4.3 Radiation 

The oscillon solution does not solve the equation of motion exactly. We have ignored terms 
of Olg'^/"^] as well as outgoing radiation. The problem of calculating the outgoing radiation 
in the small amplitude limit (not the fiat tops) has been addressed in the literature (see 
[TOl [T2]). Our intention here is to point out that for flat tops, the radiation will still be 
small. 

As shown in [21], the amplitude of the outgoing radiation can be estimated by the ampli- 
tude of the Fourier transform of the oscillon at the radiation wave number fc^ ~ \/8m (also 
see [13] )• For small amplitude oscillons, this is exponentially small ~ 6"^'''''°. Let us estimate 
what changes are expected when we move to the fiat-top oscillons. As we have seen, already 
our solutions have the form. 

Where the function 0(r, y) is independent of g. The Fourier transform of ip can be determined 
from the Fourier transform of using 

ip{t,k) = <P{T,^k). (32) 

Now (pityx) is determined entirely by a. Hence, for any given a, the Fourier transform of 
ip gets narrower as g is increased. Thus, by increasing g we can make the amplitude at the 
radiating wave number as small as we want. Note that even though the Fourier transform of 
a flat-top oscillon resembles a "sine" function, rather than a sech, this is true only for wave 
numbers near zero. Since the flat top oscillons are smooth solutions, their Fourier transforms 
still exhibit a rapid asymptotic decay. The argument in 3+1 dimensions will be similar. 
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Figure 4: In the above figure we plot N = J $^(ix for tlie oscillons in 1 + 1 (left) and 
3 + l(right) dimensions. The stability of oscillons is determined by the sign of dN/da^ where 
= g{l — uP'). Note the important difference between the curves in the two cases. While all 
oscillons are stable to long wavelength perturbations in 1 + 1 dimensions, this is not the case 
in 3 + 1 dimensions. Only those with small frequency (or equivalently, towards the flat-top 
regime) are robust. Frequency is measured in units of m. 



5 Linear Stability Analysis 



In this section we investigate whether oscillons are stable against small, localized perturba- 
tions. As discussed in the previous section, the periodic oscillon expansion, formulated in 
powers of g~^^'^, fails to solve the governing field equations and must expel radiation. In our 
stability analysis we ignore the effects of the exponentially suppressed radiation and focus 
on perturbing the oscillon profile. The main results of this section are as follows: (i) On 
the time scale of order g, 3 + 1-dimensional oscillons with large amplitudes are robust (their 
small amplitude counterparts are not) against localized perturbations with spatial variations 
comparable to the width of the oscillon. (ii) For small wavelength perturbations (compared 
to the width of of the oscillon), instabilities could exist in discrete, extremely narrow bands 
in k space. 

We now provide the details essential for reaching the above conclusions. As done pre- 
viously, we discuss the 1 -t- 1-dimensional case first, and then extend the results to 3 + 1 



dimensions. Starting with a fixed oscillon profile iposc [see Eq. (23)], we linearize about the 
oscillon by an arbitrary function x- Provided the field x remains smaller than ifosc, the 
linearized dynamics will approximately describe the perturbation. Let 



ip{t,x) = (Posc{t,x) + 6x{t,x), (33) 

where S <ti g-y^ is the amplitude of the perturbation and we keep x ~ Note that 

for a linear analysis, x must also vanish at infinity so that the perturbation 6x remains 
smaller than the original oscillon. Therefore, we restrict our analysis to spatially localized 
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perturbations. The field x then satisfies 



dtx-dix + X~^vLx + ^9^lscX = 0. (34) 

We now wish to determine if all initial conditions x remain bounded, or whether there exists 
an unstable initial profile x(0,a;)- The ip1g^{t,x),(p\g^{t^x) terms act as periodic forcing 
function^ This periodic forcing, somewhat analogous to pumping ones legs back and forth 
on a swing, may deposit energy into the field x ^i-nd consequently excite an instability. 



A complete treatment of stability may require one to solve (34) for a complete basis of 
initial conditions. Because of the spatially dependent oscillon solution, a Fourier analysis 
is difficult. With this in mind, we split the set of initial conditions into two groups. The 
first with spatial variations comparable to the size of the oscillons and another which varies 
on much shorter length scales. In the second case we can approximate the oscillon as a 
spatially constant oscillating background. This allows us to carry out a standard Floquet 
analysis. Such an analysis reveals the most dangerous instability band at /c ~ \/3, with a 
width Ak < 5f^^$Q. For large g, this becomes extremely narrow. In addition, we expect the 
time scale of these instabilities to be ~ Q^q'^- Nevertheless, one should bear in mind that 
the slow spatial variation of the oscillon could still be important. 

Now, let us look at the case where the perturbations vary on length scales comparable to 
the width of the oscillon in detail. Note that to leading order, the forcing potential: ~ ifl^^ 
is (i) (ii) smoothly varying with a natural length Xe oc y^, and (iii) oscillating with 

period 1 in the variable r = ut. The first observation implies we may use perturbation 
theory and seek an expansion for x inverse powers of g: 

X = Xo + g~^Xi + ■ ■ ■ (35) 

In the following analysis we shall work out the linear instabilities to first order in g^^ . 
However, we must keep in mind that solutions stable to order g~^, may in fact develop 
higher order instabilities over longer time scales. Since we are interested in perturbation with 
wavelengths comparable to Xe, we rescale the length x = y/gy. To capture the instability, 
however, we introduce two times: the original oscillatory time t = ut and a slow time 
T = (7"^^ The introduction of r follows from our focus on perturbations which oscillate 
near the oscillon frequency. In addition, we require a slow time T to capture variations in 



* Since (jjosc oscillates in time, the problem is essentially one of parametric resonance stability/instability. 
We are really diagonalizing the Floquet matrix - which in this case would really be an integral operator. 

^One may want to know why T — g^^t provides the important slow time scale. A back of the envelop 
calculation is as follows - consider a homogeneous background oscillating at the oscillon frequency. A naive 
perturbation series for x in powers of exhibits an oscillating term for xo, and then a term which grows 
linearly in t for xi- This means xi becomes the same order as xo when tg~^ ~ 1. Hence, there is a 
characteristic slow time T — tg~^. If we let xo be a function of both r and T, and we choose xo correctly, 
we can ensure that xi remains small for long times. In a similar fashion, if we consider higher order terms 
in the oscillon expansion, there could be additional instabilities excited over time scales g~'^, g~^ etc. 
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the perturbation. Hence, the field x = xi'^'j T, y) and the derivative dt becomes 



dt 
d1 



ujdr + g ^Ot, , 

dl + g^^[2dTdr-a'd'^] + 0[g-'] 



(36) 
(37) 



Upon substitution of Eqs. (35), (36) and (37) in Eq. (34) and collecting powers of g ^ we 
obtain 



drXo + Xo 
drXi + Xi 



0, 



-[2dTdr - a^dl - 5j - 3 cos^ r$2(?/) + 5 cos^ r$^i 
From the zeroth order equation, the most general solution for xo(t, 7", y) is 

Xo(^, T, y) = u{T, y) cos r + v{T, y) sin(r) 



Xo- 



(38) 
(39) 



(40) 



Here, u{T,y) and v(T,y) are real functions that depend on the slow time T and space y. 
Eliminating the secular terms from the right hand-side of the xi equation, we obtain 



2dTU 
2dTV 



Lv, 
-Mu, 



where the L and M are both Hermitian operators. Explicitly. 

L 

M 



-dl + a'-h\y) + l^\yl 



(41) 
(42) 



(43) 
(44) 



Since Eqs. (41) and (42) are linear, we can separate variables via u{T,y) = e^^'^u^y) 



v{T,y) 



62 



v{y): 



or equivalently 



Lv, 
-Mu, 



-LMu, 
-MLv. 



(45) 
(46) 



(47) 
(48) 
(49) 



Since both u and v are real fields and L and M are real operators, the eigenvalue^ 
must also be real. Hence, all exponents Q are either purely real or purely imaginary. Then, 
oscillon stability is guaranteed when max(f2^) < 0, or equivalently when the largest real 



^li they exist. 
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eigenvalue of —ML is negative. Determining the largest real eigenvalue of —ML can be 
done using the analysis performed by Vakhitov and Kolokolov [23] . Specifically, they exploit 
properties of the operator potentials found in L and M to show that maa;(fi^) < if and 
only if dN/da^ > 0. Here, N is the integral over all space: 

N = j ^\y)dy (50) 

and = g{l — w^). From Fig. |4| we can see that dN/da^ > for all allowed a in 1 + 1 
dimensions. Thus, 1 + 1-dimensional oscillons are stable against small perturbations with 
long wavelengths. Note that to order g~^, N = 2g^/'^Eosc in 1 + 1 dimensions (in 3 + 1 
dimensions N = 2g~^^'^Eosc-) 

The argument of [23] holds for dimensions D = 1, 2, 3. The discussion above carries over 
to 3 + 1 dimensions through the following identifications: y ^ p and dy ^ dp + {2/ p)dp 
and N = Air J ^'^{p)p'^dp. The result in 3 + 1 dimensions is in sharp contrast with that 
in 1 + 1 dimensions (see Fig. |4]). Unlike the 1 + 1-dimensional result, not all oscillons 
are robust against long wavelength perturbations. Only oscillons with large a (equivalently 
small frequency or large amplitudes) are robust. This result makes the large amplitude, 
fiat-topped oscillons in 3 + 1 dimensions particularly interesting. 

In the context of Q-balls, N is proportional to the conserved particle number and plays 
a role in the stability [20] ■ A similar interpretation might be possible here, since to leading 
order in g~^^'^, our solution is periodic in time. Finally, we note that the behavior of in 
1 + 1 and 3 + 1 dimensions can be understood heuristically. In 1 + 1 dimensions, for small a, 
the amplitude of the profile at the origin ~ a whereas the width ~ 1/a. Hence, ~ a. For 
a; — )■ etc we have increasingly wide oscillons with amplitudes ~ $c- Hence, N diverges. Now 
for 3+1 dimensions, the behavior at a — )■ is similar to the 1 + 1-dimensional case. However, 
at small a, due to the different spatial volume factor, we get ~ 1/a, therefore implying a 
nonmonotonic behavior in N. Based on a numerical analysis in the case of dilatonic scalar 
fields, it was conjectured in [12], that the stability of oscillon like configurations is related to 
the slope of the Eosc vs amplitude curve. In the large g limit, Eosc oc N and the amplitude 
oc a. Hence, their conjecture is in agreement with our analytic result. 



6 Including expansion 

In this section we consider the effects of expansion on the lifetimes and shapes of oscillons. 
We closely follow the procedure provided in [17] for the small amplitude oscillons. Here, 
applying their procedure is somewhat subtle since in the limit g ^ 1, oscillons tend to be 
very wide (oc ^/g), and the width grows without bound when a — )■ 0, Consequently, in 
these regimes it is easier to break up the oscillons due to Hubble horizon effects. Nevertheless 
we construct approximate solutions when the oscillon width is small compared to the Hubble 
horizon. 
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6.1 Including expansion in 1 + 1 dimensions 



As before, we begin with 1 + 1 dimensions and generalize to 3 + 1 dimensions. We will work 
in static deSitter co-ordinates where the metric is given by 

ds^ = -(1 - x^H^)dt^ + (1 - x^H^y^dx". (51) 

Here, if is a constant Hubble parameteij^ In these co-ordinates, the equation of motion 
becomes 

(1 - x^H^y^d^ip + 2xH^d,ip - (1 - x^H'')dl^ = -V'{ip), (52) 

where (xH) < 1. We will assume that H <^ 1 and that H = H/g where H is a small 
number. The effects of expansion can be ignored when x ^ H^^. For oscillons with widths 
satisfying Xe{a) <C H~^, the solution to the above equation is well approximated by the 
Minkowski space solution. However, in the tail of the oscillon profile we cannot ignore the 
effects of expansion. Nevertheless, taking advantage of the exponential decay of the profile 



in the tails, we can linearize Eq. (52) and obtain a solution using the WKB approximation. 



We carry out the change of space-time variables and redefinition of the field as was done 
in the nonexpanding case, Eq. dsj). Again collecting powers of g, we get 



9^01 + 01 = 0, 



(53) 



In the case of the Minkowski background, we chose an initial condition 5tv?(0,x) = 0, which 



picked out one of the two linearly independent solutions of the first equation in (53 ). However, 
in the expanding universe we need to keep the general solution 



01 (^,y) 



+ c.c. 



(54) 



where $(?/) can be complex and c.c stands for complex conjugate. The "profile" equation is 
given by 

3,_o. 5, 



{o? - [yEf] $ - (9j$ - -|$|^$ + ^l$r<I> = 



(55) 



and includes the effect of expansion through the {yHy term. We now analyze different 
regimes as seen in Fig. [s} For {yHY '^^5 the equation admits solutions identical to the 
nonexpanding case [see Eq. (23)]. In the region ye^o) <^ y ^ aH~^, where ye{o:) is the 



approximate width of the oscillon [Eq. (25)], the profile has the form 



^{y) ~ %\l ^^i^l^ exp [-ay] 



u 



(56) 



^The assumption of H being constant is for simplicity. The analysis carries over to a time dependent H 
as long as the frequency of oscillation H. 
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Figure 5: In an expanding background, if the width is small compared to H~^, the flat space 
solution is adequate for distances much less than a{y/gH)^^ . For distances larger than this, 
but still smaller than H~^, the oscillon feels the expansion, and loses energy in the form of 
outgoing waves (see inset in Fig. Distance is measured in units of the m~^. 



Since this is an exponentially decaying solution, we can ignore the nonlinear terms in the 
potential when y ^ ye{oi): 

dl^ + {{yHf-a']^^Q ye(«)«l/. (57) 

For y > aH~^, the above equation has a WKB solutioii^ in the form of an outgoing wave: 



^, , ^ 2a{l + u) 
Hy) ~ ^o\^^ — ^exp 
uHy 



AH 2 ^ 



(58) 



The amplitude of the outgoing wave was chosen using the WKB connection formula to 



1+u 



match the oscillon profile in Eq. (|56|). In terms of the original variables, we obtain 

cos(a;t) + 0[g-^^^ 
2(1 + m)« \ . 1 



1 + M cosh \2axl ^J~g\ 



where 



g^/'^H\x\u 



e cos 



ut - -Hx' 



\x\ <^a{^H) \ 
ai^Hy^ < \x\ < H-\ 



(59) 



u 



a/1 - {a/ac)' 



V~9 



(60) 



^assuming the WKB condition ff/a^ <C 1 is satisfied 
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Our solution matches that of [T7] in the hmit a <^ a^- However, as a gets larger the 
coefficient in front of the traveling wave captures the effects of the flat-top solutions. We 
will return to the above solution when we discuss the rate of energy loss by oscillons after 
considering the effects of expansion in 3 + 1 dimensions. 

We end the section by reminding ourselves of the assumptions required for this solution 
to be valid: (i) fi- > 1, (ii)il < 0{g-\ and (iii) 

Xe(a) < a{y/gH)-\ (61) 

For any H, the solution is not valid when a — )■ or ac- Also note that for a given H <^ 1, 
there always exists a g, which violates condition (iii) for all allowed a. 



6.2 Including expansion in 3 + 1 dimensions 



Now, let us include the effects of expansion for the 3 + 1-dimensional cases. The metric in 
the static deSitter co-ordinates (assuming spherical symmetry) is given by 



1 - r'H')dt' + (1 - r'Hydr' + r'dn'. 



(62) 



Following a procedure similar to the one we laid out for the 1 + 1-dimensional case, we get 
the proffie equation: 



(63) 



where p = r/y/g. The effect of expansion is included through the (pH)"^ term. For a given 
a, let the approximate width of the oscillon be Pe{c(). In the region Pe(«) <S p <^ aH~^, the 
profile has the form 



$(p) ^ /(a)^exp [-ap] 



(64) 



The lack of an analytic solution, prevents us from specifying /(a). Reverting back to the 
original variables, the solution in the spatially oscillatory regime is given by 



^{t,r) = ig-'/'aY/ 



1/2 /(«) 



e cos 



-Hr' 



a 



i^H)-' « r < H~\ (65) 



where u'^ = 1 — g ^o? . 



6.3 Energy loss due to expansion 

In this subsection, we discuss the energy loss suffered by oscillons due to the expanding 
background. As before, we start with the 1 + 1-dimensional scenario and then generalize to 
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Figure 6: The above plot shows the time it takes for an oscillon to lose most of its energy 
due to the expanding background. The plot assumes a.n H = 10~^ and g = 10. The lifetime 
is sensitive to the combination ~ g expla"^ /{gH)]. Our result is only valid when a is not too 
close to or = ^27/160, since for these value our assumption -C a{y/gH) ^ breaks 
down. Note that the lifetime should be interpreted qualitatively, since we are not allowing 
for instabilities that lead to the oscillon's abrupt disintegration. Time is measured in units 
of m~^. 



3 + 1 dimensions. The energy lost by an oscillon whose width is small compared to H ^ is 
given by 



dt 



rpx\X 



where Tj^ is the energy momentum tensor of the scalar field. We have ignored the dependence 
of the metric on x. We take X to be in the region sufficiently far away from the center. More 
explicitly, we consider X such that 



a 



iy^H)-' « |X| < H-\ 



Using Eq. (59) in Eq. (66), we get (to leading order in HX^ and g ^/^) 

dE,..r. , / /3(1 + m)3/2^ 




u 



exp[— vro; /2gH] 



(67) 



(68) 



In [T7j, a similar expression was provided in the limit where m — )■ 1. Oscillons are known 
to lose energy very slowly for a long time and then suddenly disintegrate. Our calculation 
cannot capture this disintegration. We therefore take the following as an upper bound on 
the lifetime of an oscillon (approximated by the time it take for it to lose most of its energy 
due to expansion effects) 



r<Eo 



dEosc 




dt 






(69) 



18 



where 

u = a/1 - {a/oc)'^, 

$c r, (70) 

The terms in (...) are corrections for the deviation from the sech profile. We plot the 
lifetime as a function of the in Fig. [6j We should not trust this curve in the limit 
a —J- 0, ttc, since our assumption of Xe{a) ^ {y/gH)~^ is broken here. Note the scaling with 
g: T gexp[a'^/{gH)]. 

A numerical investigation of oscillon lifetimes in a 1 + 1-dimensional expanding back- 
ground was carried out in [7]: however, no analytic calculation was provided for the lifetime. 
Our analytic results seem to be in good qualitative agreement with their paper. Note that 
Eosc ~ g~^l'^c? for small a, which tells us that the lifetime T ~ exp[i?osc/5'^'^^-f^]; is in 
agreement with the empirical formula inferred in their paper (also see Fig. 2 in [^). 

A similar calculation can be carried out in 3 + 1 dimensions. The rate of energy loss is 
then given by (to leading order in g'^^"^) 

at g 

Here, we have used Eq. (65). The lifetime is then given by 



T<-^^ge-''^^^-\ (72) 

Note that 2Eosc{ci) = g^^'^N{a) up to leading order in powers of g. We plot the approximate 
lifetime of an oscillon in a 3 + 1-dimensional scenario in Fig. [6] The lifetime is maximized at 
a < ttc- In plotting the lifetime, the numerically obtained N{a) and /(a) had errors of about 
1 percent. Again, we stress that one should think about the above curve somewhat quali- 
tatively since higher order instabilities are ignored. It would be interesting to numerically 
check the above analysis using a full 3 + 1-dimensional evolution. 



7 Conclusions 

In this paper, we pointed out what is required for scalar field potentials to support oscillons. 
We found that the spatial profiles can be very different from a Gaussian, an ansatz often made 
in the literature. In particular, we derived the nonmonotonic relationship between the height 
and the width of the oscillons, and discussed the importance of this feature for the stability 
of oscillons. We showed that the flat-topped oscillons are more stable in three dimensions 
to long wavelength perturbations as compared to their usual Gaussian counterparts. To the 
best of our knowledge, this had not been done previously in the literature in the context of 
oscillons. Oscillons could have important applications in cosmology, especially in the early 
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Universe. With this in mind, we discussed the changes in the profile, loss of energy from these 
oscillons due to expansion, and estimated their lifetimes. We provided analytic results for 
the 1 + 1-dimensional scenario, and extended analytically where possible to 3 + 1 dimensions 
and numerically otherwise. 

A number of questions related to this work require further investigation. Our expressions 
for lifetime and arguments for stability, especially in 3 + 1 dimensions in the flat-top regimes 
should be checked with a detailed numerical investigation. The question of the possible small 
wavelength, narrow band instability needs to be resolved rigorously. Recently, [21] discussed 
oscillons in the presence of gravity (an oscillaton). It would be interesting to revisit this 
problem in the context of our large energy, fiat-top oscillons. A study of oscillons emerging 
from (p)reheating-PSl[2S] like initial conditions in the early Universe is currently in progress. 
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